Ab initio coupled-cluster theory for open-shell nuclei 
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We develop a new method to describe properties of truly open-shell nuclei. This method is based 
on single-reference coupled-cluster theory and the equation-of-motion method with extensions to 
nuclei with A ± 2 nucleons outside a closed shell. We perform proof-of-principle calculations for the 
ground states of the helium isotopes 3_6 He and the first excited 2 + state in 6 He. The comparison 
with exact results from matrix diagonalization in small model spaces demonstrates the accuracy 
of the coupled-cluster methods. Three-particle-one-hole excitations of 4 He play an important role 
for the accurate description of 6 He. For the open-shell nucleus 6 He, the computational cost of the 
method is comparable with the coupled-cluster singles-and-doubles approximation while its accuracy 
is similar to coupled-cluster with singles, doubles and triples excitations. 



I. INTRODUCTION 

The nuclear shell model is the paradigm for our under- 
standing of atomic nuclei Within this model, dou- 
bly magic nuclei (with fully occupied shells for protons 
and neutrons) are particularly important because they 
are stronger bound than their neighbors and can be ap- 
proximated by a simple product of single-particle states. 
Within the nuclear shell model, doubly magic nuclei are 
the cornerstones for our understanding of entire regions 
of the nuclear chart as they can be viewed as inert cores. 
For the ab initio description of doubly magic nuclei, the 
coupled cluster (CC) method - based on particle-hole ex- 
citations of a reference Slater determinant that obey the 
linked cluster theorem - is particularly well suited and 
arguably one of the most efficient methods Simi- 
lar remarks hold for nuclei that differ from doubly-magic 
nuclei by one nucleon; such nuclei still exhibit a simple 
structure and a single Slater determinant is a good refer- 
ence state. The structure of all other nuclei is more com- 
plicated and requires the superposition of many product 
states and correspondingly large model spaces. 

For the light p-shell nuclei, various ab initio meth- 
ods [HHl yield virtually exact results for realistic Hamil- 
tonians. For heavier systems, one typically relies on ap- 
proximations. Here, coupled-cluster theory is an ideal 
compromise between accuracy on the one hand and com- 
putational cost on the other. This method has been ap - 
plied to various problems in nuclear structure @, H| fl2| — 

In this paper, we use equation-of-motion (EOM) tech- 
niques within the coupled-cluster method, abbreviated 
to EOM-CC hereafter, for the description of nuclei that 
differ from closed-shell references by two nucleons. This 
extension of the coupled-cluster method is useful for two 
reasons. First, it significantly enlarges the set of nuclei 
that can be accessed within coupled-cluster theory. For 
the oxygen isotopes 14_28 0, for instance, all nuclei ex- 
cept 19 differ from closed-subshell references by two 



neutrons or less. Similar comments apply to isotopes of 
helium and calcium. Questions related to the evolution of 
shell structure 0, Qj| could thus be addressed from first 
principles. Second, the coupled-cluster method yields a 
similarity transformed Hamiltonian (see Eq. © below) 
for a doubly-magic nucleus. The Hamiltonians for one 
and two nucleons attached to this doubly magic core pro- 
vide us with effective single-particle energies and an ef- 
fective two-body interaction, respectively. These matrix 
elements could enter the construction of effective shell- 
model interactions [2(J that are the basis for large scale 
shell- model calculations 0, |22| ■ 

Within the EOM-CC methods, the equations for ex- 
cited states and one particle attached/removed are well 
known in quantum chemistry Jj, f2a-[27| , and have also 
been applied to atomic nuclei (2 81 - 13 If . However, the cor- 
responding equations for two particle attached/removed 
have seen very few applications in quantum chem- 
istry [32l [33| and are new in nuclear physics. In this 
article, we will thus extend the range of the EOM-CC 
methods to include open-shell nuclei with A = ±2 nu- 
cleons outside a closed shell core. We present here the 
necessary formalism for deriving such equations, includ- 
ing the pertinent diagrams and algebraic equations. To 
our knowledge, these details have not been presented else- 
where. The results from our EOM-CC calculations are 
compared with full configuration interaction (FCI) calcu- 
lations for helium isotopes, demonstrating the accuracy 
of this approach. 

This paper is organized as follows. In SecHH we give a 
brief overview of the equation-of-motion method within 
coupled-cluster theory. The extension of this method to 
two valence nucleons outside a closed-shell core is pre- 
sented in Sec. Ill Bl We discuss and present our results 
in Sec. Mil Section IIVI contains our conclusions and an 
outlook for future work. 
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II. COUPLED-CLUSTER THEORY AND 
EQUATIONS-OF-MOTION FOR NUCLEI 

A. Single-reference coupled-cluster theory 

In this section we introduce the Hamiltonian that 
enters our calculation, together with a brief review of 
single reference coupled-cluster theory. We keep the 
presentation limited to those details that are required 
for the derivation of the two particle attached/removed 
(2PA/PR-EOM-CC) amplitudes presented in Sec. iHBl 
The interested reader is referred to [H for details. 

We use the intrinsic Hamiltonian 
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Here A is the number of nucleons, and iij is the nucleon- 
nucleon interaction. We will limit ourselves to two-body 
interactions only. In second quantization, the Hamilto- 
nian can be written as 



* = E 



j ^{P<l\\rs)ala\a s a r 



(2) 
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pqrs 



The term (pq\\rs) is a shorthand for the matrix elements 
(integrals) of the two-body part of the Hamiltonian of 
Eq. (JTJ, Vi1i r and s represent various single-particle 
states while e v q stands for the matrix elements of the one- 
body operator in Eq. ([TJ. Finally, operators like a q and 
a p create and annihilate a nucleon in the state q and p, 
respectively. These operators fulfill the canonical anti- 
commutation relations. 

In single-reference coupled-cluster theory, the many- 
body ground-state |\&o) is given by the exponential 
ansatz, 



|*o) =exp(f)|$ 



(3) 



Here, |<I>o} is a single- reference Slater determinant, and 
T is the cluster operator that generates correlations. 
The operator T is expanded as a linear combination of 
particle-hole excitation operators 



T = T 1 +T 2 



(4) 



where T n is the n-particle-n-hole excitation operator 

f » = E Ct"< • • ■ oL°*n • ■ • «U ■ (5) 

We use throughout this work the convention that the in- 
dices ijk . . . denote states below the Fermi level (holes), 
while the indices abc . . . denote states above the Fermi 
level (particles). For an unspecified state, the indices 



pqr . 



are used. The amplitudes t° 



will be detcr- 



singles and doubles approximation we truncate the clus- 
ter operator as 



T « Tccsd = Ti + Ti 



(6) 



which defines the coupled-cluster approach with singles 
and doubles excitations, the so-called CCSD approxima- 
tion. The unknown amplitudes result from the solution 
of the CCSD equations given by 



m\B\^) = 0, 
<*#|£|*o) = 



The term 



H = exp (-T)H exp (T) = [H exp 



(7) 



(8) 



is the similarity transform of the normal-ordered Hamil- 
tonian. The state | < i ) .fj > ") is a Slater determinant that 
differs from the reference |4>o) by holes in the orbitals 
ij . . . and by particles in the orbitals ab . . . . The sub- 
script C indicates that only connected diagrams enter. 

Once the if and tfj amplitudes have been determined 
from Eq. ([7]), the correlated ground-state energy is given 
by 



E cc = ($ |i?|$o}+£o 



(9) 



mined by solving the coupled-cluster equations. In the 



Here, Eq denotes the vacuum expectation value with 
respect to the reference state. The coupled-cluster 
equations |J7J| show that the reference state |$o) is an 
eigenstate of the similarity transformed Hamiltonian (|8]) 
within the space of lp-lh and 2p-2h excitations. 



B. Equation-of- motion coupled-cluster theory for 
open-shell nuclei 

In this work, our focus is on the development of 
coupled-cluster theory for truly open-shell nuclei, i.e., 
systems where no single reference can be constructed 
without breaking symmetries (such as rotational invari- 
ance). One could apply the CCSD method to the de- 
formed (symmetry breaking) Hartree-Fock ground state 
of an open-shell nucleus. However, the restoration of an- 
gular momentum requires more than singles and doubles 
cluster excitations (see for example Ref. [15|]) and is com- 
putationally expensive. Here, we wish to stay within the 
computationally inexpensive CCSD scheme. 

Open-shell systems can be computed with multi- 
reference methods. In such an approach, many refer- 
ence wave functions are included and treated on an equal 
footing. However, the loss of mathematical simplicity 
and transparency, and problems related to intruder states 
make these multi-reference approaches difficult to pursue. 
For a detailed discussion, we refer the reader to Ref. [13] ■ 
Equation-of-motion methods (see [26|, [2?| for recent re- 
views) avoid these problems as they exhibit the trans- 
parency and computational simplicity of single-reference 
coupled-cluster theory. 
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Within the EOM-CCSD approach, the states of the 
A ± 2 open-shell nuclei are computed from the ground 
state of the ^4-body system as 



< A) ) 



Ri A±2) 



exp(f)|$ ). (10) 



Here, i?^ 4±2 ' is a particle removal or particle addition 
operator that generates an (A ± 2)-body state from the 
A-body coupled-cluster wave function. The label /i iden- 
tifies the quantum numbers (energy, angular momentum, 
...) of the state of interest. 

The operator R^, and the energies En of the states of 
interest solve the eigenvalue problem [J, l23H27j 

(M^)c\$o)=u J »R<f ± »\$o). (11) 

Here, the expression [HR\f ±2 ^)c denotes all terms that 
connect the similarity transformed Hamiltonian H with 
the excitation operator i?^ 4±2 '. The energy difference 
Wj, = — Eq is the excitation energy of the state (x in 
the nucleus A ± 2 with respect to the ground state of the 
reference nucleus with mass A. 

The operators relevant for this work are (we drop 
the label fi for convenience) 

R (A+2) = \ E r ab aU + \ £ rfa\ afe + ... (12) 

ba iabc 



g ijk u a 



ClkCLjCLi 



,(13) 



ijka 



where the unknown amplitudes r (subscripts and super- 
scripts dropped) can be grouped into a vector that solves 
the eigenvalue problem of Eq. (jTTJ) . The operator ([12")) 
consists of a 2p-0h term, a ip-lh term, and in general 
up to an (A + 2)p-Ah term. In this work, we will trun- 
cate the operator ([12"]) at the 3p-lh level. Clearly, this 
truncation will only be a good approximation for states 
in the (^4 ± 2)-body system that have a relatively simple 
structure built on the A-body nucleus. We will intro- 
duce two different truncations for the particle attached 
and the particle removed method, and identify them by 
the number of particle-hole excitations kept in the oper- 
ator (fT2"|) . A truncation after the first term in Eq. (TT21 is 
referred to as 2PA-EOM-CCSD(2p-0/i), while the trun- 
cation after the second term is denoted as 2PA-EOM- 
CCSD(3p-lh). Similarly for 2PR-EOM-CCSD, we will 
use the abbreviations 2PR-EOM-CCSD(0p-2/i) and 2PR- 
EOM-CCSD(lp-3h) for truncations after the first and 
second term in Eq. (fT3")) . respectively. Table Q] shows the 
excitation operators used in these truncation schemes. 

We construct the matrix (i.e., the connected part of 
HR) of the eigenvalue problem of Eq. (jTTJ) diagrammat- 
ically. As usual, lines directed upwards represent parti- 
cle states, while lines directed downwards represent hole 
states [i]]- The horizontal lines represent the operators 
and we use a heavy and a wiggly line too differentiate 
the two operators in the composite diagrams. Table |TT] 
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Table I: Definition of the EOM excitation operators for the 
two particles attached (removed) method, using a truncation 
at both the 2-particle-0-hole (0-particle-2-hole) and 3-particle- 
1-hole (l-particle-3-hole) level. The operator R^-oh defines 
2PA-EOM-CCSD(2p-0h), the operator R^ lh defines 2PA- 
EOM-CCSD(3p-l/i), the operator R A ~ 2 2h defines 2PR-EOM- 
CCSD(0p-2/i) while the operator R A ~\ defines 2PR-EOM- 
CCSD(lp-3fo). These operators enter the eigenvalue problem 
of Eq. dill). 



Amplitude 



Diagram 



ijk 




Table II: Diagrams corresponding to the excitation operators 
defined in Table [I] Upward directed lines denote unoccupied 
orbitals (particle states) and downward directed lines denote 
occupied orbitals (hole states). The heavy horizontal line rep- 
resents the operator vertex, used to distinguish the R opera- 
tors from the similarity transformed Hamiltonian H. 



shows the diagrams corresponding to the r amplitudes, 
while the matrix elements of the similarity-transformed 
Hamiltonian are represented by the diagrams shown in 
Table UTT1 These elements are well known and computed 
from the corresponding contractions of the cluster opera- 
tor © with the Hamiltonian of Eq. ([1} (see, for instance, 
Ref . [J| ) . They result from the construction of the oper- 
ator ([5]) after the CCSD equations (O have been solved. 

Table IIVI shows the algebraic expressions for the ma- 
trix elements of the similarity-transformed Hamiltonian 
H. For notational efficiency, some intermediate objects 
(x) that are common among several of the matrix ele- 
ments, are defined separately in Tabic [V] In the numer- 
ical implementation, the storage of the similarity trans- 
formed Hamiltonian requires some memory. However, 
this modest cost in memory yields a significant reduction 
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Matrix element 



Diagram 



(i\H\a) 



(a\H\b) 



m\j) 



(ai\H\bc) 



(ij\H\ka) 



(ab\H\cd) 



(ij\H\kl) 



(ia\H\bj) 



(ab\H\ci) 



(ia\H\jk) 



(abc\H\dei) 



(ija\H\klm) 



{aib\H\cdj} 



(ija\H\kbl) 















Table III: Diagrams of the matrix elements of the similarity- 
transformed Hamiltonian H. The horizontal wiggly line rep- 
resents the operator vertex. The top three diagrams repre- 
sent the one-body matrix elements of H , the last four dia- 
grams represent the three-body matrix elements of H, and 
the remaining diagrams denote two-body matrix elements of 
H, respectively. The corresponding algebraic expressions are 
shown in Table HVl 



in computational cycles. For a detailed analysis we refer 
the reader to Refs. [H, [36j]. 

In a diagrammatic language, the left-hand-side of the 
eigenvalue problem (|TT|) consists of all topologically dif- 
ferent diagrams that result from connecting a diagram 
from Table HI1 with a diagram in Table [TTTT Figures Q] and 
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(ija\H\klm) 
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Table IV: Algebraic expressions for the matrix elements of 
the similarity-transformed Hamiltonian in terms of the clus- 
ter amplitudes tl and tfj, the matrix elements (ij\\ef) of 
the two-body interaction, and the one-body matrix elements 
fq — £q + ^2;{pi\ v \qi) of the normal-ordered Hamiltonian. 
The permutation operator P pq permutes the indices p and q, 
and we define P(pq) = 1 — P pq , P{pq, r) = 1 — P pr — P qr , and 
P(p, qr) = 1 — Ppq — P pr . The intermediates \ are dehned in 
Table M 



Intermediate 


Expression 


xt 


+ (a7n\\be)t^ n 


xtl 


(at || be) — ^ (7ni\\bc) 




(ij\\ka) + \ (ij\\ea)tl 


X"l°j 


(ia\\bj) + ^(ai\\eb)tj 


x'Tj 


X"^ + |<oi||e6)*5 - IC^J?' 




X'faj - ftm^T + i<^l|e6>*"TO 




(a6||cz) + ^- (a&| |ce)£| 




(to||jfc> + §<i«||e/>*^ + P{jk)t" ]X "T k - \t a m H) 



Table V: The intermediates that enter the construction of the 
similarity-transformed Hamiltonian H in Table HVl All other 
terms are defined in Table IIVI 



[2] show the diagrams of (HR)c for the two truncations 
of the 2PA-EOM-CCSD method. Let us briefly discuss 
some of these diagrams. 

For 2PA-EOM-CCSD(2p-0/i), the relevant diagrams 
correspond to the matrix element {§ ah \(HR)c\&o), i-e., 
they have two outgoing particle lines and consist of con- 
tractions of the similarity transformed Hamiltonians with 



r ab . As an example, consider the diagram of Fig. 1(b) 
It results from contracting the diagram of r ab (cf. 



Ta- 



(a) 



(b) 



(a) 



(b) 



Figure 1: Diagrams corresponding to the matrix element 
(<f> ab \{HR)c\$o) for the 2PA-EOM-CCSD(2p-0/i) amplitude 
equation. All diagrams are constructed by contracting a di- 
agram from Table (TT] with a diagram from Table IIHI Only 
diagrams that satisfy the topological form of (& ab 1 (H R)c\$o) , 
with two external particle lines in the upper part of the dia- 
gram and no external lines in the bottom part of the diagram, 
are selected. 



Figure 3: Diagrams corresponding to the matrix element 
(*y|(HJS)c|*o> for the 2PR-EOM-CCSD(0p-2/i) amplitude 
equation. All diagrams are constructed by contracting a di- 
agram from Table [II] with a diagram from Table IIHI Only 
diagrams that satisfy the topological form of ($i 3 \(HR)c\&o), 
with two external hole lines in the upper part of the diagram 
and no external lines in the bottom part of the diagram, are 
selected. 





Figure 2: Diagrams corresponding to the matrix elements 
($ ab \{HR) c \®o) (a-d) and ($f c \(HR)c\$o) (e-k) for the 
2PA-EOM-CCSD(3p-l/i) amplitude equation. 



Figure 4: Diagrams corresponding to the matrix elements 
(*y|(5A)c|*o> (a-d) and (<f>? jk \(HR) c \® ) (e-k) for the 
2PR-EOM-CCSD(lp-3h) amplitude equation. 



bleHH, with the diagram of (ab\H\cd) (cf. Table | 
For 2PA-EOM-CCSD(3p-l/i) additional diagrams of the 
form ($ ab \(HR)c\§o) enter because the amplitude rf bc 
is also permitted. In addition, diagrams corresponding 
to the matrix element (&f bc \(H R)c\$o) (i.e., diagrams 
with three outgoing particle lines and one incoming hole 
line) enter. The diagram in Fig. |2(h)] for instance, is con- 
structed by contracting the diagram element representing 
r abc j n Table mj w it n t ne diagram element representing 
(ab\H\cd) in Table HH 

Let us turn to the diagrams for two-particle removal. 
Figures [3] and 0] show the diagrams of (HR)c for the 
two truncations of the 2PR-EOM-CCSD method. Here, 
one needs topologies of the form ](HR)c\^o) and 



($° fc | (HR)c\ $o) using the r,j and r"- fe diagrams in Table 

eh 

The algebraic expressions corresponding to these dia- 
grams are derived according to the standard rules, see 
for example Ref . Q , and shown in Table I VII The com- 
putational cost is n\n Q for the 2PA-EOM-CCSD(3p-lh) 
method, and the most expensive diagram is shown in 
Fig. |2(h)| Likewise, the most expensive dia gram for the 

' and re- 
n„ is the 



4(i) 



2PR-EOM-CCSD(Dj-3/i) is shown in Fig 
quires of the order of n^n* operations. Here, 
number of unoccupied orbitals (equal to the size of the va- 
lence space) , and n a is the number of occupied orbitals in 
the reference state. For comparison, the computational 
costs of the single-reference CCSD and CCSDT meth- 
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Matrix element Expression 

2PA-EOM-CCSD(2p-0fr): 

{<S> ab \(HR) c \$o) P(ab)H b r ae +\Hf f r ef 

2PA-EOM-CCSD(3p-lfa): 

{$ ab \(HR)c\$o) P(ab)H b r ae + lS$r ef 

($f c \(HR) c \$o) P(a,bc)H b Jr^ +P(ab, c)H c e rf e 

-UTr abc +\P{ab,c)H? f r e / c 
+P(ab,c)HS c K bE +\Hfftr ef 
+±P(a,bc)H%\ c r%f 

2PR-EOM-CCSD(0p-2/i.): 

($y|(i?A)c|$o> -P{ij)H?r im +\H™ n r mn 

2PR-EOM-CCSD(lp-3fc): 

($y|(#A)c|$o> ~P{ij)H™r m + \H% n r mn 

-\-H\ e Tj i j m -^Piij^H j e Ti rnn 
($Z jk \(HR)c\$o) -P(i,jk)H^r lm +H%rl jk 

—P(ij,k)H k v i j m -\--^P(ijjk)H^j r mnfe 

* jr \ l 'J]™J £1 ek 1 tjm i 2 ijk '17171 
2 J^)^jek ^imn 



Table VI: Algebraic expressions for the 2PA/2PR-EOM- 

CCSD(2p-0/i/0p-2/i) and 2PA/2PR-EOM-CCSD(3p-l/i/lp- 
3h) approximations. All terms are denned in Table ITVl 



Method for asymmetric eigenvalue problems, see for ex- 
ample Ref. [37[ and references therein. This method is 
based on repeated applications of the matrix- vector prod- 
uct (HR)c- Specifically, our numerical implementation 
uses the ARPACK software [38] package. The expres- 
sions in Table I VII can thus be used to solve the eigen- 
value problem directly. In this paper, we employ the 
m-scheme basis for the vectors R = ( r ab ; r^ bc ) of the 
2PA. Within this scheme, we are limited to small model 
spaces. However, in this work we are mainly interested 
in testing the newly developed methods and in gauging 
the accuracy of the employed cluster truncation through 
comparisons with exact diagonalization. Note that the 
similarity-transformed Hamiltonian H exhibits the sym- 
metries of the underlying Hamiltonian. In the m-scheme 
basis, we can classify states by their projections J z of 
angular momentum and T z of isospin and their parity. 
Although the solutions will have a good total angular 
momentum J 2 , we will not be able to exploit this sym- 
metry in the m-scheme basis. 

The FCI method we employ [3!| is limited to rela- 
tively small model spaces. Note, however, that we are 
also working on an angular-momentum-coupled imple- 
mentation gg of 2PA/2PR-EOM-CCM. This will allow 
us to exploit rotational symmetries and give us access to 
much larger model spaces. These spaces are well beyond 
the reach of present full diagonalization methods. 



III. RESULTS 



ods are n\r? and n^n 3 , respectively. This shows that 
the 2PA / 2PR-EOM-CC methods developed in this work 
are relatively inexpensive from a computational point of 
view, and comparable to the cost of CCSD. Note that the 
three-body matrix elements of the similarity-transformed 
Hamiltonian are not stored. Instead, they are calculated 
when needed and do not contribute significantly to the 
overall cost of the calculations. For a detailed analysis 
we refer the reader to Refs. [35l [36j . 

The eigenvalue problem of Eq. (fTTj) for the two-particle 
attached system thus becomes 



($ ab \(HR) c \1> Q ) = 
($f c \(HR)c\$o) = 



(14) 



Here, the left-hand-side is a linear function of the "vec- 
tor" R = ( r ab ,rf c ) of amplitudes (see Table EJ) and 
constitutes a matrix- vector product. Note that the two 
equations are coupled and constitute a single eigenvalue 
problem for the 3p-lh truncation. Likewise, we find for 
the two-particle removed problem 



(^ 3 \(HR)c\^o) 
($? jk \(HR) c \*o) 



ijk 



(15) 



We are usually only interested in the few lowest eigen- 
values of Eq. dm) . For this purpose, we use the Arnoldi 



For the proof-of-principle study we consider the he- 
lium isotopes 3 He to 6 He. Here, 3 He and 5 He are viewed 
as one neutron removed from and attached to 4 He, re- 
spectively. In a shell-model picture, 6 He is a truly open- 
shell nucleus with two valence neutrons in the P3/2 shell. 
Thus, a single-reference Slater determinant may not be 
a good starting point for coupled-cluster calculations of 
this nucleus, and it seems advantageous to describe this 
six-nucleon system as two (halo) neutrons added to the 
4 He core. We do not present results for 2 He using the 
2PR-EOM-CCSD approach; an analysis of this method 
will be presented elsewhere (ioj . 

For our calculations we use a realistic nucleon-nucleon 
(N N) p otential derived from chiral effective field theory 
[41M45 1 at order N 3 LO using interaction matrix elements 
from Ref. 45[. The matrix elements of the bare interac- 
tion employ a cutoff at A = 500 MeV. The short-range 
parts of the interaction are removed via the similarity 
renormalization group transformation [461 ] with a cutoff 
at 1.9 fin -1 . We use the spherical harmonic oscillator 
with the oscillator frequency haj = 24 MeV as our single- 
particle basis. Our model space consists of five major 
oscillator shells, with maximum orbital angular momen- 
tum Z max = 2. This results in a total of 76 single-particle 
states for neutrons and protons. We neglect three-body 
and four-body interactions. This leads to missing many- 
body physics but is not relevant for our proof-of-principle 
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computation and comparison with FCI calcualtions. 

We calculate the ground state energies of the A= 3-6 
helium isotopes. For 6 He we also compute the first 2+ 
state and the expectation values of the total angular mo- 
mentum. In addition, we also discuss the first 1 _ excited 
state, as an example where the Zp-lh truncation fails. 
We compare the equation-of-motion (EOM) approach to 
the FCI method and to three different single-reference 
coupled-cluster approximations. Recall that within the 
m-scheme, one can also compute open-shell nuclei di- 
rectly with the coupled-cluster method without resort- 
ing to EOM techniques. While such a direct approach 
may suffer from the lack of a good reference state (as 
is the case for the ground state of 6 He), the inclusion 
of more and more clusters must converge to the FCI re- 
sults. Comparing the EOM-CC approach to these single- 
reference coupled-cluster calculations allows us to gauge 
the efficiency of the various coupled-cluster approxima- 
tions. 

The direct coupled-cluster calculations employ the 
CCSD approximation described above, the CCSDT ap- 
proximation (that includes triples clusters) and the 
CCSDT- 1 approach which includes some of the 3p-3h 
clusters of the full CCSDT approximation. For the EOM 
calculations, we use the CCSD wave function of 4 He 
as the reference wave function and employ the intrin- 
sic Hamiltonian of Eq. (JTJ with A = 3,5,6 for 3 ' 5,6 He, 
respectively. We compare all results to an exact calcu- 
lation done with the FCI approach, using the same in- 
teraction and model space. The ground-state energies of 
3 " 5 He are shown in TablcEE!] For 3 He, CCSDT becomes 
an exact method and agrees with FCI. Here, the single- 
reference coupled-cluster calculations are superior to the 
EOM-CCSD approach. Evidently, the weakly bound nu- 
cleus 3 He is not well described as a neutron removed from 
the tightly bound 4 He. This means in turn that corre- 
lations beyond one-particle-two-hole excitations play a 
non-ncglig ible role. For 4 He, the EOM-CCSD approach 
is identical to CCSD. Here, triples corrections and full 
triples represent a significant improvement over CCSD, 
bringing the results close to the FCI ones. For 5 He, 
the EOM-CCSD is superior to a single-reference CCSD 
calculation and competes well with the computationally 
more expensive triples correction CCSDT-I. Clearly, the 
valence neutron in 5 He is weakly correlated with the 
strongly bound 4 He core, and the PA-EOM-CCSD ap- 
proach captures this state very well. 

We turn to the truly open-shell nucleus 6 He and 
show our results in Table 171111 Here, the 2PA-EOM- 
CCSD (3p-lh) approach is clearly superior to the single- 
reference coupled-cluster approaches, as it reproduces the 
energy and the spin of the ground state to a very good 
approximation. For the computation of the spin within 
the single-reference approaches, we compute the expecta- 
tion value (J 2 ) within the Hellmann-Feynman theorem, 
and define (J) from the relation (J 2 ) = (J) ((J) + 1). A 
direct computation within the CCSD or CCSDT approxi- 
mations can only be based on a symmetry-breaking refer- 



^Hc ^Hc 5 He 

CCSD -6.624 -27.468 -22.997 

CCSDT-1 -6.829 -27.600 -23.381 

CCSDT -6.911 -27.619 -23.474 

EOM-CCSD -6.357 -27.468 -23.382 

FCI -6.911 -27.640 -23.640 



Table VII: Ground-state energies (in MeV) for 3 He, 4 He 
and 5 He, calculated with coupled-cluster methods trun- 
cated at the 2-particle-2-hole (CCSD) level, 3-particle-3-hole 
(CCSDT) and a hybrid (CCSDT-1) where a small subset of 
the leading diagrams in CCSDT are included. For the EOM- 
CCSD approach, truncations has been made at the 1-particle- 
2-hole level, the 2-particle-2-hole level, and the 2-particle-l- 
hole level for 3 He, 4 He and 5 He respectively. The energies are 
compared to the exact full configuration interaction (FCI). 





o| 


2+ 


0+ {J) 


2 + 


CCSD 


-22.732 


-20.905 


0.78 


2 


CCSDT-1 


-24.617 


-21.586 


0.25 


2 


CCSDT 


-24.530 


-21.786 


0.01 


2 


2PA-EOM-CCSD(2p-0/i) 


-21.185 


-18.996 





2 


2PA-EOM-CCSD(3p-l/i) 


-24.543 


-21.634 





2 


FCI 


-24.853 


-21.994 


{) 


2 



Table VIII: Energies (in MeV) for the ground state and first 
excited state of 6 He and the expectation value of the total 
angular momentum, calculated with coupled-cluster methods 
truncated at the 2-particle-2-hole (CCSD) level, 3-particle-3- 
hole (CCSDT) and a hybrid (CCSDT-1) where the 3-particle- 
3-hole amplitudes are treated perturbatively. The 2PA-EOM- 
CCSD results are calculated with a truncation at the 2- 
particle-O-hole (2PA-EOM-CCSD(2p-0/i) level and at the 3- 
particle-l-hole (2PA-EOM-CCSD(3p-l/i) level. All energies 
are compared to full configuration interaction (FCI) results. 



ence state. Clearly, 2-particle-2-hole excitations (CCSD) 
cannot restore the symmetry, and we find that (J) = 0.78 
for the ground state of 6 He. Adding 3-particle-3-hole 
excitations (CCSDT approximation) almost restores the 
rotational symmetry, but some correlation energy is still 
missing. In the 2PA-EOM-CCSD approach, however, the 
rotational symmetry is preserved throughout the calcu- 
lation, and we obtain a very good approximation of the 
energy at a relatively low computational cost. As ex- 
pected, the 2PA-EOM-CCSD (2p-0h) approach is less 
accurate than the 3p-lh approximation, since it is unable 
to account simultaneously for the correlations within the 
three-body system consisting of the two valence neutrons 
and the 4 He core. The ground state and the first excited 
2 + state of 6 He are both dominated by a configuration 
with two neutrons in the p 3 / 2 orbit. This is consistent 
with the shell-model picture of this nucleus. For the ex- 
cited 2+ state, the CCSD, CCSDT and CCSDT-1 meth- 
ods result in the correct value of the angular momentum 
due to the choice of reference state. 

Let us define the fraction <t coit of the correlation energy 
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E< 



2 PA 



En 



Efci — E 



(16) 



Here E 2 p A is the energy from 2PA-EOM-CCSD, E FC i 
is the energy from the FCI, while Eq is the energy ex- 
pectation of the uncorrelated reference state |$o)- Both 
E2PA and Epci are shown in Table IVIII1 and Eq — 
— 16.807 MeV. We also compute the norm 



^ = £l 



■ib\2 



£ 



ib\2 



(17) 



ab 



abi 



The normalized squared weights 



ab 

p\ ee Af- 1 £ |r? 

abi 



b 1 2 



61 2 



(18) 



fulfill Pi + p| — 1 an d measure the importance of the 
2p-0h and 3p-lh amplitudes, respectively. 

Table IIXI shows the relative correlation energies <7 corr 
and the relative weights pi and p 2 of 2PA-EOM- 
CCSD(3p-l/i) for the three lowest states with quantum 
numbers + , 1~ and 2 + of 6 He, respectively. We see 
that 2PA-EOM-CCSD accounts for more than 90% of 
the correlation energy for the + and 2 + states, and 
that most of the weight is carried by the 2p-0h ampli- 
tude. The 1~ state, however, is not very accurately 
reproduced and much of the correlation energy is lack- 
ing. For this state, the 2PA-EOM-CCSD(3p-l/i) energy 
is Ei- = —20.95 MeV and deviates considerably from 
the full CI result E x - = -23.26 MeV. Consistent with 
this picture is the large weight carried by the 3p-lh am- 
plitudes. For a converged computation, one would pre- 
sumably also need to include 4p-2h or higher clusters. 
Inspection shows that the J v = 1~ state is dominated 
by two neutrons in the 0p 3 / 2 single-particle state, but 
with an additional 1-particle-l-hole excitation of either 
a proton or a neutron to the OP3/2 state or the 0pi/ 2 or- 
bit. These 3p-lh configurations are energetically favored 
compared with a configuration with one neutron in the 
OP3/2 state and one in the 0^5/2 state, a configuration 
which can also give a 1" state. This explains why this 
state is dominated by the rf bc amplitudes. Note finally 
that we cannot expect a separation of the center-of-mass 
motion from the intrinsic dynamics in the small model 
space we considered (UH^- Thus, the low- lying 1~ state 
might also exhibit considerable admixtures with spurious 
center-of-mass excitation. 

The results of this section show that states in open- 
shell nuclei that exhibit a simple structure imposed on 
a correlated core can be well described by the EOM 
method. Moreover, the EOM wave function preserves 
the symmetries of the Hamiltonian and is expected to 
be useful in the computation of matrix elements besides 
the energy. Taking into account the low computational 





Ccorr 


Pi 


P2 




0.96 


0.84 


0.16 


ir 


0.64 


0.34 


0.66 


o + 

'I 


0.93 


0.81 


0.19 



Table IX: The relative correlation energy <T CO rr defined in 
Eq. (|16[) for the lowest states with quantum numbers J w — + , 
1~ and 2 + in 6 He, respectively. The relative weights pi and p2 
of the 2p-0h and 3p-lh amplitudes, respectively, are defined 
in Eqs. JT8j>. 



cost of the EOM-CC methods as compared to the full 
CCSDT approach, our EOM approach is clearly well 
suited for these selected states. For states where more 
complicated particle-hole excitations are prominent, the 
various truncation schemes discussed here are insufficient 
and we will need additional correlations in the EOM op- 
erator to reach satisfactory results. The ground state of 
3 He and the first excited 1~ state of 6 He discussed above, 
are examples in case. 



IV. CONCLUSIONS AND FUTURE 
PERSPECTIVES. 

We have developed and implemented the equation- 
of-motion coupled-cluster method and performed micro- 
scopic calculations of helium isotopes with up to two 
valence nucleons outside the closed-shell alpha parti- 
cle. The comparison with full configuration interaction 
calculations shows that the equation-of-motion coupled- 
cluster methods yield accurate results for sufficiently sim- 
ple states. The open-shell nucleus 6 He, for instance, 
can be viewed and computed as two weakly correlated 
neutrons attached to the correlated core of 4 He. The 
proof-of-principle calculations were performed in a re- 
duced model space for a comparison with results from 
exact diagonalizations. We are working on implement- 
ing our formalism for the two-particle attached and re- 
moved equation-of-motion coupled-cluster methods in an 
angular momentum coupled basis. This will allow us to 
employ much larger model spaces. With this improve- 
ment, the first-principles computation of semi- magic nu- 
clei, single-particle energies, and effective two-particle in- 
teractions for the nuclear shell-model can be addressed. 
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